Quantum Lifshitz point in the infinite dimensional Hubbard model 
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We show that the Gutzwiller variational wave function is surprisingly accurate for the compu- 
tation of magnetic phase boundaries in the infinite dimensional Hubbard model. This allows us 
to substantially extend known phase diagrams. For both the half-hypercubic and the hypercubic 
lattice a large part of the phase diagram is occupied by an incommensurate phase, intermediate 
between the ferromagnetic and the paramagnetic phase. In case of the hypercubic lattice the three 
phases join at a new quantum Lifshitz point at which the order parameter is critical and the stiffness 
vanishes. 
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The Hubbard model was originally introduced to study 
ferromagnetism in strongly correlated metals [J 0, 
This phenomenon as prototypically realized in Fe, Co 
and Ni is one of the oldest phenomena investigated in 
solid state theory. A related problem is incommensu- 
rate spin-density-wave order as manifested in Cr and its 
alloys Also various transition metal oxides show in- 
commensurate magnetic phases which are often accom- 
panied by charge order, like cuprates nickelates @ 
and manganites[7|. 

Despite decades of investigation not much is known 
about the magnetic phase diagram of the Hubbard 
model. The simplest treatments [l, Q partition the zero 
temperature phase diagram in three regions: antiferro- 
magnetic (AFM) close to n = 1 particles per atom, fer- 
romagnetic (FM) at large interaction and far from n — 1 
and paramagnetic (PM) at small interaction and/or close 
to n = 0, 2. More sophisticated treatments include in ad- 
dition incommensurate (IC) phases I(J 11 1. 

Infinite dimensional lattices offer an unique opportu- 
nity to study the competition between PM and FM keep- 
ing the problem tractable and, at the same time, retain- 
ing much of the physics expected in three dimensional 
lattices [l2j]. Here we investigate the Hubbard model 
in the limit of infinite dimension T> using the Gutzwiller 
variational wave function (GWF) [13]. Instabilities of the 
FM and PM ground state are systematically studied as a 
function of momentum, doping, and interaction strength 
usin g a random-phase-approximation (RPA) like expan- 
sion 
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principle the model can be solved ex- 
actly in this limit by using dynamical mean-field theory 
(DMFT) which maps the problem to an impurity prob- 
lem amenable of numerical solution (l2| . Limitations on 
the numerical algorithm, however, have prevented for an 
extensive zero temperature investigation of the phase di- 
agram. A study by Unrig on the stability of the FM 
phase in an infinite dimensional generalization of the fee 
lattice, the so called half-hypercubic (hhc) lattice, is one 
of the few cases where the T = self-consistent DMFT 



problem has been solved exactly [141 ]. 

Here we show, comparing with exact results when 
available, that the celebrated GWF is surprisingly accu- 
rate for the determination of magnetic phase boundaries 
in infinite dimensional lattices. In addition we signifi- 
cantly extend the computation of the T = magnetic 
phase diagram to regions in parameter space yet poorly 
explored by DMFT methods. More specifically, for the 
hhc, we show that the GWF FM instability line agrees 
with the exact computation by Uhrig 14j . In addition we 
present the first systematic investigation of the stability 
of the PM state in this lattice and find that the PM and 
FM never meet but are separated by an IC phase (Fig. [5]). 
A similar systematic study of the hypercubic lattice (he) 
shows also a large region of, yet poorly studied, IC or- 
der (Fig. [3]). In contrast to the hhc case, PM and FM 
phases have a common phase boundary. The IC phase, 
the FM phase and the PM phase join in a quantum ver- 
sion of the multicritical Lifshitz point (LP). Classical LP 
are interesting because they are associated with unusual 
critical exponents [HI • We anticipate unusual critical be- 
havior for a quantum LP (QLP) shall it occur at physical 
dimensions. 

We start from the Hubbard Hamiltonian [H, 0, El 



H — tijC ia Cj C 



U 



(1) 



where c[ a (cj CT ) destroys (creates) an electron with spin 
a at site i. U is the on-site Hubbard interaction and tij 
denotes the hopping parameter between sites i and j. 

The energy of the model is evaluated within the GWF 
for polarized and unpolarized phases. In the limit T> — > 
oo the Gutzwiller approximation (GA) to the variational 
problem becomes exact |20j which greatly simplifies the 
computations. The stability of these solutions is studied 
by an RPA analysis at different momenta. By construc- 
tion, the limits of stability thus obtained correspond to 
the point where a more stable variational solution would 
be found in a fully unrestricted computation of the GWF 
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FIG. 1: (color online) Magnon dispersions for the (fully po- 
larized) ferromagnet in the half-hypercubic lattice for S = 0.3. 
At this doping the FM solution is stable for U > 0.98Ubr- 



energy. This requires a rotationally invariant formula- 
tion of the GA energy functional [2l| and to take into 
account the change in the double occupancy for small 
quasistatic changes of the charge and spin distribution 
as in Vollhardt's GA-based Fermi liquid computation 
[22I ]. Basically we compute the dynamic spin suscepti- 
bility x 9 M = —fa Je lult (TS+(t)SZ q (0)) in both para- 
and FM states which can be obtained from a density ex- 
pansion of the Gutzwiller energy functional. Details of 
the formalism are given in Ref. [171 ]. For simplicity we 



neglect macroscopic phase separation [9], which, in 
any case, would be frustrated if the long-range Coulomb 
interaction were taken into account [l9j |. 

We consider the hypercubic lattice and the half- 
hypercubic lattice. The latter can be obtained from the 
hypercubic lattice by removing all the even sites. For the 
hhc case we consider hopping restricted to Uj = t/T> for 
nearest neighbor sites and Uj = t/(2T>) for next-nearest 
neighbor sites For the hypercubic lattice we keep 
only nearest-neighbor hopping tij = tj \J2T>. In the lat- 
ter case it has been shown [23| that all the momentum 
dependence of correlation functions is contained in the 
factor 



1 V 



(2) 



which can be parametrized by a value between 1 and —1 
corresponding to a scan along the zone diagonal from q = 
(0, 0) to q = (tt, 7r), respectively. On the other hand 
for the half-hypercubic lattice the Brillouin zone is cut in 
half and the momentum dependence is contained in the 
factor ?7q which varies from 1 to corresponding to a scan 
from q = (0, 0) to q = (tt/2, 7r/2), respectively [53]. 

The hhc lattice is non-bipartite, therefore antiferro- 
magnetism is frustrated. In addition, it has a density of 
states with a square root divergence at a lower band edge 
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FIG. 2: (color online) Full line: limit of stability of the FM 
phase within the Gutzwiller wave function for the hhc lattice. 
The asterisks are the exact result from Ref. Il4l . In the inset a 
scaled representation of the FM stability limit is shown with 
U r ed = U/(U + Ubr) and Ubr = 6.86/;. Dashed line: limit of 
stability of the PM phase. The arrows indicate the value of 
?7q parameterizing the momentum of the unstable mode. 



that makes ferromagnetism stable for large U and large 
8. Here 5 = 1 — n is the density deviation from half- 
filling and n is the density. For convenience we define 
the Brinkmann-Rice critical U value Ubr = — 8e with e 
the noninteracting kinetic energy at half- filling [22| . 

Fig. [T] shows the transverse magnon dispersion com- 
puted on top of a saturated FM solution at 6 = 0.3 in 
the hhc for different values of U which as expected is 
linear (quadratic) in 1 — ?7 q (q) for small q and shows a 
Goldstone mode at q = 0. As U is decreased the exci- 
tation energy becomes negative first at q = and then 
also at finite momentum signaling an instability due to a 
change of sign in the stiffness of the magnetization. Since 
this is a transverse excitation the ferromagnet tends to 
become a spiral. 

In Fig. [2] we show the limit of stability of the saturated 
FM phase thus obtained (full line) and the exact com- 
putation by Uhrig (asterisks). Only in the region where 
bound states appear for small 6 [141 ] significant devia- 
tions appear as shown in the inset which indicates that 
the GWF energy is quite accurate. 

We also show in Fig. [5] the stability limit of the PM 
phase (dashed line). The first excited state of the para- 
magnet is a triplet with z-component of the spin to = 
— 1,0, 1. Spin rotational invariance requires that the lon- 
gitudinal (to = 0) and transverse (to = — 1, 1) excita- 
tions are degenerate and indeed all three become soft at 
the boundary at some finite q. In the figure we indicate 
by arrows the value of 7/q for the unstable modes. The 
instability occurs at 7y q — * for small 8, consistent with 
the wave vector q = (%/2, 7r/2) and gradually shifts 
to rjq = 1 for 8 = 1 as expected for the instability toward 
the FM phase (q = 0). 

Both instability lines indicate that there is an IC phase 
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FIG. 3: (color online) Limit of stability of the FM (full line) 
and PM (dashed line) phases for the he lattice. Arrows indi- 
cate the value of rjq for the unstable mode of the PM phase. 
The stars correspond to the limit of stability of the FM phase 
obtained in Ref. within an approximate solution of the 
DMFT equations. Dotted lines are levels of constant magne- 
tization. Energies are in units of Ubr = 6.38i. 



intermediate between the FM phase and the PM phase. 
This does not show up in the DMFT numerical study of 
Ref. |24| because it is restricted to q = instabilities. 
The present FM boundary practically coincides with the 
T = extrapolated results of Ref. 1241 while a substantial 
portion of the PM range of Ref. [24| is instead occupied 
by IC phases in the present computation. 

Encouraged by these results which show the accuracy 
of the Gutzwiller variational energy to determine the 
phase boundaries, we consider the he lattice. The limit of 
stability of the FM phase and of the PM phase are shown 
in Fig. [31 Since the he lattice has perfect nesting the in- 
stability towards antiferromagnetism (rjq = —1) sets in 
at U/t = for 5 = @, [Tlj]. For 6 > an IC phase 
emerges with ?? q gradually increasing. Decreasing U at 
finite <5 the FM phase becomes unstable at the full line 
due to a change of sign in the stiffness of the magnetiza- 
tion as for the hhc lattice. Finally at the QLP the PM 
instability mode becomes uniform (7y q = 1) and one has 
a transition to a weak (partially polarized) FM phase. 

The dotted lines in Fig. [3] are level curves of constant 
magnetization m. The density of states of the T> = oo hy- 
percubic lattice is a Gaussian. Because of the Gaussian 
tails an infinite splitting of the minority and majority 
bands, i.e. U — > oo, is required to produce a fully polar- 
ized state. For finite U one has a magnetization of the 
FM state that increases from zero as S decreases from the 
PM boundary (dashed line) and becomes logarithmically 
saturated for small S. 

The limit of stability of the FM phase in Fig. [3] is in 
good agreement with the FM phase boundary for the 
hypercubic lattice obtained in a DMFT study based on 
a large U expansion and the no crossing approximation 
as impurity solver (asterisks) 25]. This again points to a 



good accuracy of the GWF. 

The magnetic phase diagram of the T> = oo he lattice 
has been studied at finite temperature by Freericks and 
Jarrell using DMFT (2(|. The finite temperature insta- 
bilities are towards a commensurate AFM and become 
IC only at the lowest temperature which does not allow 
for a reliable T = extrapolation of the PM-IC phase 
boundary. If we restrict to the AFM instability of the 
GWF we find an instability line in the U — S plane (not 
shown) which is in good agreement with the T — ex- 
trapolated phase diagram of Ref. [26| . This is probably 
due to the fact that the extrapolation is dominated by 
the high temperature commensurate data. In the present 
phase diagram the IC phase occupies a larger region. 

At the QLP the IC, PM and FM phases meet. This 
behavior is the quantum analogue of the classical LP ob- 
served in MnP where the same phases meet but at finite 
temperature and in the presence of an external field [13] • 

On general grounds regarding the proximity of T> = 3 
systems to infinite dimensional systems [12J we expect the 
phase diagram of Fig. [3] to resemble that of the cubic 
lattice. It is then worth speculating on the behavior of a 
QLP shall it occur at finite T>. In analogy with a classical 
LP also for a QLP we expect anomalous critical behavior. 

The T = PM-FM transition and PM-IC transition 
are the first studied examples of quantum critical behav- 
ior [II SHIM El III. The transition is characterized 

by a dynamical exponent z — 3 (z = 2) for the PM-FM 
(PM-IC) and the upper critical spatial dimension is 4 — z. 
Therefore in many cases of interest the theory is Gaus- 
sian. It was pointed out that in the case of FM transi- 
tions because the magnetization is a conserved quantity 
long-range interactions appear and the critical exponents 
change [3l|. 

At the QLP the critical behavior is expected to be 
different. Neglecting long-range interactions the effective 
action close to the QLP can be written in analogy to a 
classical LP and the FM QCP [29|: 



S 



d v q 

r + cq 2 + Dq A + 7 ^ ) 0(q, w„).0'(-q, -W n ) 



(3) 



d v x I dT[0(x,r).0(x,r)] 2 . 



Here </>(x, r) is a vector order parameter defined in space 
and imaginary time, ^>(q, u> n ) is the Fourier transform 
and uj n is a bosonic Matsubara frequency. D, u and 7 
are positive constants. The parameters r and c are linear 
functions of <5 and U close to the quantum critical point. 
At mean-field level the QLP is defined by c = 0, r = 0. 
At the FM-PM boundary r = and c > whereas at 
the PM-IC phase boundary c < and r = c 2 /(4D) and 
the unstable mode has momentum q = y— c/ (2D). The 
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FM-IC boundary is characterized by c = 0, r < thus 
c changes sign at the boundary mimicking the change of 
sign in the stiffness found in RPA. 

As mentioned above, far from the QLP the transitions 
have a dynamical exponent z — 2,3. At the QLP the dy- 
namical exponent changes to z = 5 and the upper critical 
dimension becomes 8 — z — 3. Therefore, if a QLP exists 
at low T> the quartic term is not irrelevant and the the- 
ory becomes non-Gaussian or Gaussian with logarithmic 
corrections in marked contrast with Hertz theory. A full 
analysis must take into account singular contributions to 
the coefficients and goes beyond our present scope. 

One can also envisage an AFM QLP at which AFM, 
PM and IC phases meet. In this case the damping term 
becomes momentum independent implying z = 4 with 
an upper critical dimension T> = 4. Cuprates may be 
not so far from this phenomenology since at low temper- 
atures and as a function of 5 they show a transition from 
an AFM state to a paramagnet/superconductor which in 
addition often shows IC behavior [5] . Even more relevant 
for our discussion may be the low temperature behavior 
of overdoped cuprates where the IC states turns into a 
paramagnet. It has been argued[34] that the system is 
also close to a FM instability which suggests the proxim- 
ity to a IC-PM-FM QLP. Another interesting system is 
CeCug-ajAua; for which an AFM QLP where only one di- 
rection becomes soft has been proposed [35| . In all these 
cases, however, disorder effects appear to be quite rele- 
vant. Another system where one may hope to find a QLP 
are ultracold atoms in optical lattices. In this fascinating 
realization of the Hubbard model one can fine tune the 
parameters in the Hamiltonian which may allow for a full 
exploration of the phase diagram [36| ■ 

To conclude, our comparison with DMFT studies 
shows that the GWF, which played and still plays a fun- 
damental role in understanding strong correlation, is a 
very convenient and accurate tool to study the T> — > oo 
magnetic phase diagram of the Hubbard model. Wc 
find that quite generically incommensurate phases ap- 
pear as intermediate phases between strongly polarized 
FM phases and PM phases. This is interesting because, 
as mentioned in the introduction, IC phases are observed 
in a variety of systems in physical dimensions. For the 
he lattice we find that the weakly polarized FM and the 
PM phase have a common phase boundary and both meet 
with the IC phase at a new QLP at which the magnetic 
order parameter is critical and the stiffness simultane- 
ously vanishes. We speculate that an analogous QLP 
may exist in lower dimensional systems and we antici- 
pate anomalous quantum critical behavior. 
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